A CMV-Based Eigensolver for Companion Matrices
نویسندگان
چکیده
In this paper we present a novel matrix method for polynomial rootfinding. We approximate the roots by computing the eigenvalues of a permuted version of the companion matrix associated with the polynomial. This form, referred to as a lower staircase form of the companion matrix in the literature, has a block upper Hessenberg shape with possibly nonsquare subdiagonal blocks. It is shown that this form is well suited to the application of the QR eigenvalue algorithm. In particular, each matrix generated under this iteration is block upper Hessenberg and, moreover, all its submatrices located in a specified upper triangular portion are of rank two at most, with entries represented by means of four given vectors. By exploiting these properties we design a fast and computationally simple structured QR iteration which computes the eigenvalues of a companion matrix of size n in lower staircase form using O(n2) flops and O(n) memory storage. So far, this iteration is theoretically faster than the fastest variant of the QR iteration for companion matrices in customary Hessenberg form. Numerical experiments show efficiency and accuracy of the proposed approach.
منابع مشابه
Lightweight 4x4 MDS Matrices for Hardware-Oriented Cryptographic Primitives
Linear diffusion layer is an important part of lightweight block ciphers and hash functions. This paper presents an efficient class of lightweight 4x4 MDS matrices such that the implementation cost of them and their corresponding inverses are equal. The main target of the paper is hardware oriented cryptographic primitives and the implementation cost is measured in terms of the required number ...
متن کاملCompression of unitary rank-structured matrices to CMV-like shape with an application to polynomial rootfinding
This paper is concerned with the reduction of a unitary matrix U to CMV-like shape. A Lanczos–type algorithm is presented which carries out the reduction by computing the block tridiagonal form of the Hermitian part of U , i.e., of the matrix U +UH . By elaborating on the Lanczos approach we also propose an alternative algorithm using elementary matrices which is numerically stable. If U is ran...
متن کاملToward High Performance Divide and Conquer Eigensolver for Dense Symmetric Matrices
This paper presents a high performance eigensolver for dense symmetric matrices on multicore architectures. Based on the well-known divide and conquer (D&C) methodology introduced by Cuppen, this algorithm computes all the eigenvalues of the symmetric matrix. The general D&C can be expressed in three stages: (1) Partitioning into subproblems, (2) Computing the solution of the subproblems and (3...
متن کاملDarboux transformations for CMV matrices
We develop a theory of Darboux transformations for CMV matrices, canonical representations of the unitary operators. In perfect analogy with their self-adjoint version – the Darboux transformations of Jacobi matrices – they are equivalent to Laurent polynomial modifications of the underlying measures. We address other questions which emphasize the similarities between Darboux transformations fo...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- SIAM J. Matrix Analysis Applications
دوره 36 شماره
صفحات -
تاریخ انتشار 2015